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ABSTRACT 

The quasar 3C 279 was the target of an extensive multiwavelength monitor- 
ing campaign from January through April 2006. An optical-IR-radio monitoring 
campaign by the Whole Earth Blazar Telescope (WEBT) collaboration was or- 
ganized around Target of Opportunity X-ray and soft 7-ray observations with 
Chandra and INTEGRAL in mid- January 2006, with additional X-ray coverage 
by RXTE and Swift XRT. In this paper we focus on the results of the WEBT 
campaign. 

The source exhibited substantial variability of optical flux and spectral shape, 
with a characteristic time scale of a few days. The variability patterns throughout 
the optical BVRI bands were very closely correlated with each other, while there 
was no obvious correlation between the optical and radio variability. After the 
ToO trigger, the optical flux underwent a remarkably clean quasi-exponential 
decay by about one magnitude, with a decay time scale of ~ 12.8 d. 
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In intriguing contrast to other (in particular, BL Lac type) blazars, we find 
a lag of shorter-wavelength behind longer-wavelength variability throughout the 
RVB wavelength ranges, with a time delay increasing with increasing frequency. 
Spectral hardening during flares appears delayed with respect to a rising optical 
flux. This, in combination with the very steep IR-optical continuum spectral 
index of a Q ~ 1.5-2.0, may indicate a highly oblique magnetic field configuration 
near the base of the jet, leading to inefficient particle acceleration and a very steep 
electron injection spectrum. 

An alternative explanation through a slow (time scale of several days) accel- 
eration mechanism would require an unusually low magnetic field of B < 0.2 G, 
about an order of magnitude lower than inferred from previous analyses of simul- 
taneous SEDs of 3C 279 and other FSRQs with similar properties. 

Subject headings: galaxies: active — Quasars: individual (3C 279) — gamma- 
rays: theory — radiation mechanisms: non-thermal 



1. Introduction 

Flat-spectrum radio quasars (FSRQs) and BL Lac objects are active galactic nuclei 
(AGNs) commonly unified in the class of blazars. They exhibit some of the most violent 
high-energy phenomena observed in AGNs to date. Their spectral energy distributions 
(SEDs) are characterized by non-thermal continuum spectra with a broad low-frequency 
component in the radio - UV or X-ray frequency range and a high-frequency component 
from X-rays to 7-rays. Their electromagnetic radiation exhibits a high degree of linear 
polarization in the optical and radio bands and rapid variability at all wavelengths. Radio 
interferometric observations often reveal radio jets with individual components exhibiting 
apparent superluminal motion. At least episodically, a significant portion of the bolometric 
flux is emitted in > 100 MeV 7-rays. 46 blazars have been detected and identified with high 
confidence in high energy (> 100 MeV) 7-rays by the Energetic Gamma-Ray Experiment 
Telescope (EGRET) instrument on board the Compton Gamma-Ray Observatory (CGRO, 



Hartman et al.l 1 19991 : iMattox. Hartman. fc Reimerl 1200 ll ) . 



In the framework of relativistic jet models, the low-frequency (radio - optical/UV) 
emission from blazars is interpreted as synchrotron emission from nonthermal electrons in a 



1 For questions regarding the availability of the data from the WEBT campaign presented in this paper, 
please contact the WEBT President Massimo Villata at villataOoato . inaf . it 
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relativistic jet. The high-frequency (X-ray - 7-ray) emission could either be produced via 
Compton upscattering of low frequency radiation by the same electron s responsible for the 



synchrotron emission (leptonic jet models; for a recent review see, e.g., iBottcherl l2007al ) . or 
due to hadronic processes initiated by relativistic p rotons co-accelerated with the electrons 
(hadronic models, for a recent discussion see, e.g. 
20031 ). 



Miicke fc Protheroe 2001; Miicke et al. 



The quasar 3C279 {z = 0.538) is one of the best-observed flat spectrum radio quasars, 
not at last because of its prominent 7-ray flare shortly after the launch of CGRO in 1991. 
It has been persistently detected by EGRET each time it was observed, even in its very 
low quiescent states, e.g., in the winter of 1992 - 1993, and is known to vary in 7-ray 
flux by roughly two orders of magnitude (IMaraschi et al.l 1 1994 ; IWehrle et al.l Il998l ) . It has 
been monitored intensively at radio, optical, and more recently also X-ray frequencies, and 



has been the subject of intensive mult iwavelength campaigns (e.g.. IMaraschi et al.l 11994 
Hartman et al1ll996l : IWehrle et aJll998h . 



Also at optical wavelengths, 3C 279 has exhibited substantial variability over up to two 
orders of magnitude (R ~ 12.5-17.5). Variability has been observed on a variety of different 
time scales, from years, down to intra-day time scales. The mo st extreme variability pattern s 
include intraday variability with flux decays of < 0.1 mag /hr (IKartaltepe fc Balonekl 120071 ) . 
Observations with the International Ultraviolet Explorer in the very low activity state of 
the source in December 1992 - January 1993 revealed the existence of a thermal emission 
component, possibly related to an accretion disk, with a lumi nosity of Ljju ~ 2 x 10 46 erg s~ ] 



i f this component is assumed to be emitting isotropically (jPian et al.l Il999l ). iPian et al. 



( 119991 ) have also identified an X-ray spectral variability trend in archival ROSAT data, 
indicating a lag of ~ 2 - 3 days of the soft X-ray spectral hardening behind a flux increase. 
Weak evidenc e for spectral variability was also found within the EGRET (MeV - GeV) 
energy range (INandikotkur et al.l 120071 ). At low 7-ray flux levels, an increasing flux seems 
to be accompanied by a spectral softening, while at high flux levels, no consistent trend was 
apparent. 

The quasar 3C 279 was the first object in which supe rluminal motion was discovered 
( Whitney et al. 1971 ; Cotton et al. 1979 ; Unwin et al. 19891). Characteristic apparent speeds 
of individual radio comp onents range up to /3 app ~ 17 (jCotton et al.l Il979l ; iHoman et al. 
20031 ; IJorstad et al.ll2004l ). indicating pattern flow speeds with bulk Lorentz factors of up to 
r ~ 17. Radio jet components have occasionally been observed not to follow straight, ballistic 



„ 

Homan et al. 


2003; 


Jorstad et al. 


2004) 
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Ojha et al.ll2004j ; iLister fc Homanll2005l ; iHelmboldt et al.ll2007l ). indicating that the magnetic 
field might be predominantly perpendicular to the jet on those length scales. 

A complete compilation and modeling of a ll available SEDs simul taneous with the 11 
EGRET observing epochs has been presented in lHartman et al.l (12001a[). The modeling was 



done u s ing the time-dependent lep tonic (SSC + EIC) model of lBottcher. Mause. fc Schlickeiser 



Bednarek 


1998: 


Sikora et al. 


2001 



Moderski et al.ll2003l ) concluding that the X-ray - soft 7-ray portion of the SED might be 
dominated by SSC emission, while the EGRET emission might require an additional, most 
likely external-Compton, component. The resulting best-fit parameters were consistent with 
an increasing bulk Lorentz factor, but decreasing Lorentz factors of the ultrarelativistic elec- 
tron distribution in the co-moving frame of the emission region during 7-ray high states, 
as compared to lower 7-ray states (IHartman et al.ll2001al ). However, such an interpretation 
also required changes of the overall density of electrons, and the spectral index of the in- 
jected electron power-law distribution, which did not show any consistent trend with 7-ray 
luminosity. 



Hartman et al.l (l2001bl ) have investigated cross correlations between different wavelength 
ranges, in particular, between optical, X-ray, and 7-ray variability. In that work, a general 
picture of a positive correlation between optical, X-ray and 7-ray activity emerged, but no 
consistent trends of time lags between the different wavelength ranges were found. 

The discussion above illustrates that, in spite of the intensive past observational efforts, 
the physics driving the broadband spectral variability pr operties of 3C 279 are still rather 
poorly understood. For this reason, ICollmar et al.l (l2007bl ) proposed an intensive multiwave- 
length campaign in an optical high state of 3C 279, in order to investigate its correlated 
radio IR optical - X-ray - soft 7-ray variability. The campaign was triggered on Jan. 
5, 2006, when the source exceeded an R-band flux corresponding to R = 14.5. It involved 
intensive radio, ne ar-IR (JHK), and optical monitoring by the Whole Earth Blazar Telescope 



(WEBlE 



see, e.g. 



Raiteri et al. 200 



op 



Villata et al.l 120071 . and references therein) collabora- 
tion through April of 2006, focusing on a core period of Jan. and Feb. 2006. In order to 
illustrate the source's behaviour leading up to the trigger in January 2006, previously un- 
published radio and optical data from late 2005 are also included in the analysis presented in 
this paper. X-ray and soft 7 observations were carried out by all instruments on board the 
International Gamma-Ray Astrophysics Laboratory (INTEGRAL) during the period of Jan. 
13 - 20, 2006. Additional, simultaneous X-ray coverage was obtained by Chandra and Swift 



• fottp : //www . to . astro . it/blazars/webt 
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XRT. These observations were supplemented by extended X-ray monitoring with the Rossi 
X-Ray Timing Explorer (RXTE). In this paper, we present details of the data collection, 
analysis, and results of the WEBT (radio IR optic al) campaign. Prelimi nary results of 
the mu ltiwavelength campaign have been presented in IColfmar et al.l (l2007af ) and IBottcher 
(l2007al lbL and a final, comprehensive report on the result of the entire multiwavelength 



campaign will appear in ICollmar et al.l (l2007bl ) . 



Throughout this paper, we refer to a as the energy spectral index, F v [Jy] oc v~ a . A 
cosmology with Q m = 0.3, = 0.7, and H = 70 km s _1 Mpc -1 is used. In this cosmology, 
and using the redshift of z — 0.538, the luminosity distance of 3C 279 is dz, = 3.1 Gpc. 



2. Observations, data reduction, and light curves 

3C 279 was observed in a coordinated multiwavelength campaign at radio, near-IR, 
optical (by the WEBT collaboration), X-ray (Chandra, Swift, RXTE PC A, INTEGRAL 
JEM-X), and soft 7-ray (INTEGRAL) energies. The overall timeline of the campaign, along 
with the measured long-term light curves at radio and optical frequencies is illustrated in 
Fig. HJ Simultaneous X-ray coverage with all X-ray / soft 7-ray telescopes mentioned above 
was obtained in the time frame Jan. 13 - 20, as indicated by the gray shad ed area in Fig. 



[ H Det ailed results of those high-energy observations will be presented in ICollmar et al. 



J2007bl ) 



Table [T] lists all participating observatories which contributed data to the WEBT 
campaign. In total, 25 ground-based radio, infrared, and optical telescopes in 12 countries 
on 4 continents contributed 2173 data points. 



2.1. Optical and near-infrared observations 



The observing strategy and data analysis followed to a large extent the standard pro- 
cedure for the optical data reduction for WEBT campaigns which is briefly outlined below. 
For more informa t ion on standard data redu c tion procedu r es for WEBT campaigns se e 
also: 



Villata et al.l (12 000) 



Raiteri et a 



Villata et all fl2004al ibl); iRaiteri et al.l (120051 ) 



(200 



Bottcher et al.l (120051 ) 



Villata et al l (120021) ; IBottcher et all (120031 ): 



It had been suggested that, optimally, observers perform photometric observations al- 
ternately in the B and R bands, and include complete (U)BVRI sequences at the beginning 
and the end of each observing run. Exposure times should be chosen to obtain an opti- 
mal compromise between high precision (instrumental errors less than ~ 0.03 mag for small 
telescopes and ~ 0.01 mag for larger ones) and high time resolution. If this precision re- 




720 730 740 750 760 770 780 790 800 810 820 830 840 

JD - 2453000 

Fig. 1. — Timeline of the broadband campaign on 3C 279 in 2006, including the optical and 
radio light curves during the entire campaign period. The gray shaded area indicates the 
period of the INTEGRAL, Chandra, and Swift observations. 
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Table I. List of observatories that contributed data to the WEBT campaign 



Observatory 


Specifications 


frequency / filters 




Radio 


Metsahovi, Finland 


14 m 


37 GHz 


70 


Medicina, Italy 


32 m 


5, 8, 22 GHz 


32 


Noto, Italy 


32 m 


8, 22 GHz 


6 


RATAN-600, Russia 


576 m (ring) 


1, 2.3, 5, 8, 11, 22 GHz 


138 


Crimean Astr. Obs., Ukraine (RT-22) 


22 m 


36 GHz 


7 


Infrared 


Roque (NOT), Canary Islands 


2.56 m 


J, H, K 


3 


Optical 


Abastumani, Georgia (FSU) 


70 cm 


R 


127 


ARIES, Naintal, India 




R 


63 


Bclogradchik, Bulgaria 


60 cm 


V, R, I 


75 


BOOTES-1, Spain 


30 cm 


R 


151 


Catania, Italy 


91 cm 


U, B, V 


33 


Crimean Astr. Obs., Ukraine 


70 cm 


B, V, R, I 


47 


SMARTS, CTIO, Chile 


90 cm 


B, V, R 


33 


Kitt Peak (MDM), Arizona, USA 


130 cm 


U, B, V, R, I 


190 


Kitt Peak (MDM), Arizona, USA 


240 cm 


R 


77 


Michael Adrian Obs., Germany 


120 cm 


R 


9 


Mt. Lcmmon, Arizona, USA 


100 cm 


B, V, R, I 


214 


Mt. Maidanak (AZT-22), Uzbekistan 


150 cm 


B, V, R, I 


44 


Osaka Kyoiku, Japan 


51 cm 


V, R, I 


494 


Roque (KVA), Canary Islands 


35 cm 


R 


75 


Roque (NOT), Canary Islands 


256 cm 


U, B, V, R, I 


7 


SARA, Arizona, USA 


90 cm 


B, V, R, I 


242 


Tenagra, Arizona, USA 


81 cm 


B, V, R, I 


19 


Torino, Italy 


105 cm 


B, V, R 


3 


Tuorla, Finland 


103 cm 


R 


84 
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quirement leads to gaps of 15 - 20 minutes in each light curve, we suggested to carry out 
observations in the R band only. Observers were asked to perform bias and dark correc- 
tions as well as flat-fielding on their frames, and obtain instrumental magnitudes, applying 
either aperture photometry (using IRAF or CCDPHOT) or Gaussian fitting for the source 
3C 279 and four recommended comparison stars. This calibration has then been used to 
convert instrumental to standard photometric magnitudes for each data set. In the next 
step, unreliable data points (with large error bars at times when higher-quality data points 
were available) were discarded. Our data did not provide evidence for significant variability 
on sub-hour time scales. Consequently, error bars on individual data sets could be further 
reduced by re-binning on time scales of typically 15-20 min. The data resulting at this 
stage of the analysis are displayed in the top panel of Fig. [TJ 

In order to provide information on the intrinsic broadband spectral shape (and, in 
particular, a reliable extraction of B - R color indices) , the d ata were then de-reddened using 



the Galactic Extinction coefficients of iSchlegel et al.l (119981 ). based on As = 0.123 mag and 
E(B - V) = 0.029 mag3 

Possible contaminations of the optical color information could generally also arise from 
contributions from the host galaxy and the optical - UV emission from an accretion disk 
around the central supermassive black hole in 3C 279. However, these contributions are not 
expected to be significant in the case of our campaign data: Assuming absolute magnitudes 
of M v ~ —23 and Mp ~ — 21 for typical quasar host galaxies at z ~ 0.5 (e.g. iFloyd et al. 



2004 ; IZakamska et al.ll2006l ). their contribution in the V and B band, respectively, at the 
distance of 3C 279 would be V g& \ ~ 19.5 and B ga \ ~ 21.6, respectively. These are at least 
about four magnitudes fainter than the actually measured total B and V magnitudes during 
our campaign, and thus negligible. The possible contrib ution of an accret ion disk can be 



estimated on the basis of the thermal component for which iPian et al.l (119991 ) found evidence 
in IUE observations during the 1992/1993 low state of 3C 279. Their best fit to this compo- 
nent suggests U ~ 18.6 and B ~ 21 (and much fainter contributions at lower frequencies), 
which corresponds to a contribution of < 2.5 % to the total B and U magnitudes measured 
during our campaign. Therefore, both the host galaxy and the accretion disk contribution 
are neglected in our further analysis. 

The only infrared observations obtained for this campaign were one sequence of JHK 
exposures taken on January 15, 2006, with the 2.56 m NOT on Roque de los Muchachos on 
the Canary Island of La Palma. The resulting fluxes are included in the SED displayed in 
Fig. El 



http : //nedwww. ipac . caltech . edu/ 
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2.1.1. Optical light curves 

The optical (and radio) light curves from December 2005 to April 2006 are displayed in 
Fig. [TJ The densest coverage was obtained in the R band, and the figure clearly indicates 
that the variability in B, V, and I bands closely tracks the R-band behaviour. The coverage 
in the U band was extremely sparse and does not allow any assessment of the U-band light 
curve during our campaign. Therefore, the U-band will be ignored in the following discussion, 
and we will describe the main features of the variability behavior based on the R-band light 
curve. 

The optical light curves show variability with magnitude changes of typically < 0.5 mag on 
time scales of a few days. The most notable exception to this relatively moderate variability 
is the major dip of the brightnesses in all optical bands right around our coordinated X-ray / 
soft 7-ray observations in January 2006 (~ JD 2453742 - 2453770). In the R-band, the light 
curve followed an unusually clean exponential decay over 1.1 mags, in 13 days, i.e., a slope 
of dR/dt = 0.085 mag/day or a flux decay as F(t) = F(t ) e~ <yt ~ t °^ Td with a decay constant 
of Td = 12.8 d. Only moderate intraday deviations on a characteristic scale of < 0.1 mag/d 
are superposed on this smooth exponential decay. 

In contrast to the smooth decline of the optical brightness during January 6 - 20, 2006, 
the subsequent re-brightening to levels comparable to those before the dip, appears much 
more erratic and involves a remarkably fast rise by ~ 0.5 mag within ~ 1 d on Jan. 27 (JD 
2453763). Unfortunately, the detailed shape of this fast rise was not well sampled in our 
data set. If this was a quasi-exponential rise with a slope of dR/dt ~ 0.5 mag/d, it would 
correspond to a rise time scale of r r ~ 2.2 d. 



2.2. Radio observations 

At radio frequencies, 3C 279 was monitored using the 14 m Metsahovi Radio Telescope 
of the Helsinki University of Technology, at 37 GHz, the 32-m radio telescope of the Medicina 
Radio Observatory near Bologna, Italy, at 5, 8, and 22 GHz, the 32-m antenna of the Noto 
Radio Observatory on Sicily, Italy, at 8 and 22 GHz, the 576-m ring telescope (RATAN-600) 
of the Russian Academy of Sciences, at 1, 2.3, 5, 8, 11, and 22 GHz, and the 22-m RT-22 
dish at the Crimean Astrophysical Observatory, Ukraine, at 36 GHz. 



The Metsahovi data have been reduced with the standard procedure described in lTerasranta et al 



( 119981 ) . The resulting 37 GHz light curve is reasonably well sampled during the period mid- 
January - mid-March 2006. Inspection by eye in comparison to the optical light curves 
displayed in Fig. [1] appears to indicate that the optical and 37 GHz light curves follow simi- 
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lar variability patterns with a radio lead before the optical variability by ~ 5 days. However, 
a discrete cross-correlation analysis between the R-band and 37 GHz radio light curves did 
not reveal a significant signal to confirm this suggestion. 

Also included in Fig. [I] are the radio light curves at 5, 8, and 22 GHz. For details of 
the analysis of data from the Medicina and Noto radio observatories at those frequencies, 



sec 



Bach et al.l ( 120071 ). As already apparent in Fig. [fl most of the data at frequencies below 
37 GHz were not well sampled on the < 3-months time scale of the 2006 campaign, and any 
evidence for variability did not show a discernable correlation with the variability at higher 
(radio and optical) frequencies. 

All radio data contributed to this campaign are stored in the WEBT archive (see 
http : / /www . to . astro . it/blazar s/webt/ for information regarding availability of the data) . 
Radio data at all observed frequencies have been included in the quasi- simultaneous SED 
for Jan 15, 2006, shown in Fig. [8j Given the generally very moderate radio variability at 
frequencies below 37 GHz, a linear interpolation between the two available data points near- 
est in time to Jan. 15, 2006, was used to construct an estimate of the actual radio fluxes at 
that time. 



3. Optical spectral variability 

In this section, we will describe spectral variability phenomena, i.e. the variability of 
spectral (and color) indices and their correlations with monochromatic source fluxes. We will 
concentrate here on the optical spectral variability as indicated by a change of the optical 
color. In particular, our observing strategy was optimized to obtain a good sampling of 
the B R color index as a function of time. Since our data did generally not indicate 
substantial flux changes on sub-hour time scales, we extracted B R color indices wherever 
both magnitudes were available within 20 minutes of each other. Fig. [2] shows the B R 
color history over the entire campaign, compared to the B and R band light curves. Overall, 
there is no obvious correlation between the light curves and the color behavior of the source 
on long time scales. However, two short-term sequences attracted our attention: There 
is a sequence of brightness decline accompanied by a spectral hardening (declining B R 
index) around JD 2453750 (Jan. 14, 2006), and another sequence of a brightness increase 
accompanied by a spectral softening around JD 2453790 - 2453806 (Feb. 23 - Mar. 11). 
However, we caution that incomplete sampling, in particular in the B-band may introduce 
spurious effects. 

A similar trend was recently observed in a multiwavelength WEBT campaign on the 
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Fig. 2. — Light curve of the B and R magnitudes and the B R color index of 3C 279 over 
the duration of the entire campaign. 



-14- 



I 



4e12 h 

f--.. 
3e12 h 



N 

X 

~3 



2e12 - 



1e12 - 



R V B U 

1 1 1 r- 



♦ Dec. 23 
b— a Jan. 15 

* Jan. 24 
v v Feb. 20 

Mar. 25 



1 1- 



•3 

•IL 



4e14 



5e14 6e14 7e14 8e14 
v[Hz] 



Fig. 3. — Snap-shot optical (BVRI) continuum spectra for 5 epochs during our campaign. 
All measurements for each individual spectrum were taken within < 20 min. of each other. 
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quasar 3C 454.3 (jVillata et al.ll2006l ). where it could be interpreted as a "little blue bump 
due to the unresolved contribution from optical emission lines in the ~ 2000 - 4000 A wave- 



length range in the rest frame of the quasar, in particular from Fe II and Mg II (jRaiteri et al. 



20071 ). In order to test whether such an interpretation would also be viable in the case of 



3C 279, we have compiled several simultaneous snap-shot optical continuum (BVRI) spectra 
at various brightness levels (Fig. [3]). All measurements for each individual spectrum dis- 
played in Fig. [3] have been taken within < 20 min. of each other. Only the spectrum of 
Jan. 15 shows a significant deviation from a pure power-law in the B band; and in that case, 
there was simultaneous U band coverage, which matched a straight power-law extrapolation 
of the VRIJHK spectrum. If the spectral upturn towards the blue end of the spectrum were 
due to an unresolved Mg II / Fe II line contribution, it should emerge even more clearly in 
the Jan. 24 spectrum, which is characterized by a lower optical continuum flux level than 
the Jan. 15 spectrum. Therefore, the compilation of spectra in Fig. [3] does not provide any 
support for the existence of an essentially non-variable continuum component at the blue 
end of the spectrum. For this reason, we believe that the color changes that we found in our 
data are in fact intrinsic to the blazar jet emission. 

Fig. H] illustrates our impression from Fig. [21 that there is no clear overall trend 
of source (R-band) brightness with optical spectral hardness. However, the figure clearly 
illustrates that the scatter of the B R color index is significantly larger at larger source 
brightness, indicating that spectral variability is more likely to occur when the source is 
bright. Specifically, for R > 14.5, the data is consistent with a roughly constant value of B - 
R = 1.09, corresponding to a spectral index of a Q = 1.9 for a power-law continuum spectrum 
with i^ )0 [Jy] oc v~ a °. At brightness levels R < 14.5, spectral variability by A(B — R) < 0.35, 
corresponding to Aa Q < 0.85, is observed. 

In Fig. 0, we focus on two ~ 1-month optical flares, around mid- Jan. - mid-Feb. 
2006, and late Feb. - late March 2006, as indicated by the gray shaded segments 1 and 2, 
respectively, in Fig. [2j The number labels in the color- magnitude diagrams in Fig. [5] indicate 
the time ordering of the points. While there is no obvious trend discernible in Segment 1, 
Segment 2 suggests the presence of a spectral hysteresis pattern: The spectral softening 
around JD 2453790 - 2453806 (Feb. 23 - Mar. 11), already mentioned above, precedes the 
main brightness increase. Subsequently, the optical continuum hardens while the source is 
still in a bright state. We need to caution that due to the poor sampling of the B band 
light curve during segment 2, the significance of the tentative hysteresis found here may be 
questionable. However, such hysteresis would naturally lead to a B-band time lag behind the 
R-band, for which we do find a 3.9 a evidence from a discrete correlation function analysis 
as described in the following section. 
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Fig. 4. — Color- magnitude diagram for the entire campaign period. 
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Fig. 5. — Color-magnitude diagrams for the two time segments marked in Fig. [2j with time 
ordering indicated by the numbers in the two panels. 
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To our knowledge, such a spectral hysteresis has never been observed at optical wave- 
lengths for any flat-spectrum radio quasar. It is reminiscent of the spectral hysteresis 
occasionally seen at X-ray energies in high-frequency peaked BL Lac objects (HBLs, e.g. 



Takahashi et al.lll996t iFossati et al.ll2000l ; iKataoka et al.ll2000l ). However, the spectral hys- 



teresis observed in the X-rays of HBLs is generally clockwise (i.e., spectral hardening pre- 
cedes flux rise; softening precedes flux decline), and can be interpreted as the synchrotron 
signature of fast acceleration of ultra relativistic electrons, followed by a gradual decline on 
the radiative cooling time scale fe.g.. IKataoka et al.l |2000| ; iKusunose. Takahara. fc Lil |2000| ; 
Li fc Kusunos3 l200o[ liottcher fc ChianeJ l20Q2h . In our case, the direction of the spectral 
hysteresis is counterclockwise (i.e., spectral softening precedes the flux rise; spectral harden- 
ing precedes flux decline). Possible physical implications of such hysteresis phenomena will 
be discussed in SjHl 



4. Inter-band cross-correlations and time lags 

The result of an occasional counterclockwise hysteresis in 3C 279, as found in the pre- 
vious section, immediately suggests the existence of a characteristic time lag of higher- 
frequency behind lower-frequency va riability. In order to inv estigate this, we evaluated the 



discrete correlation function (DCF, lEdelson fc Kroliki Il988l ) between the R band and the 



other optical light curves. In our notation, a positive value of the time lag At would indi- 
cate a lag of the R-band light curve behind the comparison light curve. As mentioned earlier, 
the radio (and near-IR) light curves are too sparsely sampled and yielded no significant fea- 
tures in the DCF. Fig. [6] shows the DCF between the R band and the V band (top panel) 
and the B band (bottom panel), using a sampling time scale of At = 1 d. We have done 
the same analysis using various other values of At, which yielded results consistent with the 
ones described below. We chose to show the results for At = 1 d because they provided the 
best compromise between dense time scale sampling and reduction of error bars. 

The DCFs reveal clear correlations between the different optical wavebands, with peak 
values around 1. This confirms our previous impression from inspection by eye, that the 
variability patterns in all optical wavebands track each other very closely. 

The resulting DCFs have then been fitted with a symmetric Gaussian, F e~ ( - T ~ T °^ : /( 2 °" 2 ). 
This analysis yields non-zero offsets of the best-fit maxima, r at the ~ 3 - 4cr level. Specif- 
ically, we find r = (-1.14 ± 0.48) d for the V band and r = (-3.75 ± 0.96) d for the B 
band. This indicates a hard time lag of higher-frequency variability behind the variability at 
lower frequencies in the B-V-R frequency range. However, this trend does not continue into 
the I band. This is illustrated in Fig. where we plot the best-fit time lags as a function of 
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Fig. 7. — Best-fit time lag of the R band vs. B, V, R, and I bands, as a function of frequency. 
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photon frequency. However, we need to add a note of caution: The rather sparse sampling of 
the B- and V-band light curves leads to large error bars on the DCFs. Clearly, alternative, 
more complex representations, e.g., multiple Gaussians and/or asymmetric functions, will 
also provide acceptable fits to the observed DCFs and may lead to different quantitative 
results concerning the involved time lags. Future observations with denser B- and V-band 
sampling are needed in order to test the robustness of the trend found here. 

The hard lag found in this analysis may be physically related to the ~ 2 - 3 d time 
lag of the soft X-ray spectral in dex behind the flux in ROSAT observations of 3C 279 in 
December 1992 - January 1993 (IPian et al.lll999l ). Possible physical interpretations will be 
discussed in 361 



5. Broad-band spectral energy distribution 



The most complete broadband coverage of 3C 279 during our campaign was obtained 
on January 15, 2006. On that day, the only near-infrared (JHK) exposures were taken with 
the NOT, and simultaneous radio observations at 5, 8, 22, and 37 GHz were available as 
well. This is also within the time window of the X-ray and soft 7-ray observations, although 
in most cases, a meaningful extraction of spectral information requir ed the integration ove r 
most of the high-energy observing period, January 13 - 20, 2006 (jCollmar et al.l l2007al ). 
The total snap-shot SED composed of all available radio, near-IR, optical, X-ray, and soft 
7-ray observations around January 15, 2006 is displayed in Fig. [HJ where all near-IR and 
optical data are taken within ±1/2 hr of UT 05:00 on January 15. The figure compares the 
January 2006 SED to previous SEDs from the bright flare during the first EGRET observing 
epoc h in June 1991, the lo w state of December 1992 / January 1993 (both SEDs adapted 
from lHartman et al.ll2001al ). and a previous multiwavelen gth campaign around INTEGRAL 
AO-1 observations of 3C 279 in a low state in June 2003 (jCollmar et al.l 120041 ) . 



While the optical (presumably synchrotron) emission component clearly indicates that 
the source was in an elevated state compared to previous low states, the simultaneous X-ray 
- soft-7-ray spectrum is perfectly consistent with the low-state spectra of 1992/1993 and 
2003. This is a very remarkable result and will be discusse d in detail in a comp anion paper 
about the results of the entire multiwavelength campaign (jCollmar et al.ll2007bl ). 



The NIR - optical continuum can be very well represented by a single power-law with 
a spectral index of a Q = 1.64 ± 0.04, corresponding to an underlying non-thermal electron 
distribution with a spectral index o f p = 4.28 ± 0.0 8 , if the optical continuum is synchrotron 
emission. As already pointed out in lHartman et al.l (j2001al ) for a comparison between various 
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Fig. 8. — Simultaneous snap-shot spectral energy distributions of 3C 279 at various epochs. 
Data pertaining to the multiwavelength campaign in January 2006 are plotted with red 
triangles. All IR and optical data were taken within ±1/2 hr of UT 05:00 on Jan. 15, 2006; 
the X-ray and soft 7-ray data represent the time-averaged spectra throughout the respective 
observing windows in Jan. 2006, which included Jan. 15 for all instruments. Model fits to 
the EGRET PI and P2 SEDs are calculated us ing a time-dependent leptonic ( SSC + EC) jet 



model , and taken from lHartman et al.l (j2001al ). The June 2003 SED is from ICollmar et al. 
(l2004h . 
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EGRET observing epochs over ~ 10 years and confirmed by our campaign data on time scales 
of weeks - months (see £|3]), the optical spectral index in 3C 279 does not show any systematic 
correlation with the brightness state of the source, and Fig. |H] indicates that in our high- 
state observations of 2006, the continuum spectral slope is not significantly different from 
the slopes observed in the low states of 1992/1993 and 2003. 



6. Discussion 

In this section, we discuss some general physical implications and constraints that our 
results can place on source parameters. In the following discussion, we will parameterize 
the magnetic field in units of Gauss, i.e., B = 1 B G Gauss, and the Doppler boosting factor 
D = (r[l — j3r cos^obs]) -1 in units of 10, i.e., D = 10Di where T is the bulk Lorentz 
factor of the emitting region, /3r c is the corresponding speed, and ^obs is the observing 
angle. The characteristic variability time scale is of the order of 1 - a few days, so we write 
^°ar = l^vard- m the same sense, we parameterize the time lag of the B band behind the 
R band as r^f = 1 r^ d d with r^ d ~ 3. The observed variability time scale yields an 
estimate of the size of the emitting region, R = 10 15 i?is cm through R < c.D/(l + z) t°^. 
We find R 15 < 17/^0 

The (co-moving) energies of electrons emitting synchrotron radiation at their character- 
istic peak frequencies in the R and B bands are 



lR =3.7x 10 3 (BgD,)^ 2 

lB =4.5xl0 3 (5 GJ D 1 )- 1/2 (1) 

The steep underlying electron spectrum with p = 4.5, inferred from the steep optical 
continuum, might indicate that the entire optical spectrum is produced by electrons in the 
fast-cooling regime. This implies that the radiative cooling time scale of electrons emitting 
synchrotron radiation in the optical regime is shorter than the characteristic escape time 
scale of those electrons. The respective time scales in the co-moving frame, tg SC and t' cool can 
be written as 



C =V^ < 5.7 x10 s 77 D x tt> 

4oi =i ~7.7xl0 7 J B G 1 7 - 1 (l + fc)^ 1 s (2) 

where 77 > 1 is the escape time scale parameter (as defined in the first line of eq. [2]), and 
k is a correction factor accounting for radiative cooling via Compton losses in the Thomson 
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regime in a radiation field with an energy density u' ra d = ku' sy . Requiring that the cooling 
time scale is shorter than the escape time scale, at least for electrons emitting in the R band, 
leads to a lower limit on the magnetic field: 



B > B c = 1.3 x 10- 3 (1 + ^-V 1 (Ctd)- 2 G (3) 

which does not seem to pose a severe constraint, given the values of B ~ a few G typically 
found for other FSRQ s and also for 3C 279 from previous SED modeling analyses (e.g., 



Hartman et al.ll2001al ) 



Another estimate of the co-moving magnetic field can be found by assuming that the 
dominant portion of the time-averaged synchrotron spectrum is emitted by a power-law 
spectrum of electrons with N e (j) = n Vb1~ v for 7i < 7 < 72; here, Vb is the co-moving 
blob volume, and we use p = 4.5 as a representative value inferred from the optical continuum 
slope. The normalization constant uq —n e (l—p)/ (j 2 ~ P ~~ 7i~ P ) is related to the magnetic 
field through an equipartition parameter e# = u's/u'e (in the co- moving frame). Note that 
this equipartition parameter only refers to the energy density of the electrons, not accounting 
for a (possibly dominant) energy content of a hadronic matter component in the jet. Under 



these assumptions, the magnetic field can be estimated as described, e.g., in iBottcher et al. 



(120031 ) . Taking the vF v peak synchrotron flux / e sy at the dimensionless synchrotron peak 



photon energy e sy = E p ^ sy /(m e c 2 ) ~ 3 x 10 7 as ~ 10 10 ergs cm 2 s 1 , we find 



B> B CB ^ 1.86 D^efittUr 1 G. 



(4) 



This constitutes a more useful and realistic magnetic-field estimate than eq. [31 If, in- 
deed, the optical emission is synchrotron emission from a fast-cooling electron distribution, 
then electrons have been primarily accelerated to a power-law distribution with an injection 
index of q = p — 1 = 3.5. This is much steeper than the canoni cal spectral index of q ~ 2 .2 



-2.3 



1999 



bund for acceleration o n relativistic, parallel shocks (e.g.. iGallant. Achterberg fc Kirk 



Achterberg et a 



Qstrowski Bednarz 



20021 : 



2001 ), and could indicate an ob lique magnetic-field orientation (e.g., 



Niemiec fc Ostrowskill2004l ). which would yield a consistent pic- 



ture w ith the predominan t ly perpendicular magnetic-field orienta tion observed on parsec- 
scales (IJorstad et al.ll2004l ; lOjha et al.l I2004J : iLister fc Homanl 120051 ) . The observed hard lag 
(B vs. R) may then be the consequence of a gradual spectral hardening of the electron accel- 
eration (injection) spectrum throughout the propagation of a relativistic shock front along 
the jet. Such a gradual hardening of the electron acceleration spectrum could be the con- 
sequence of the gradual bui ld-up of hydromagnetic turbulence through the relativistic two- 
stream instability (see, e.g.. ISchlickeiser et al.ll2002l ). This turbulence would harden the rel- 
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atiyistic electron distribution via 2nd-order Fermi acceleration processes (jVirtanen fc Vainio 



20051 ). Such a scenario would imply a length scale for the build-up of turbulence of Ar 



ct^ DT/(1 + z) ~ 5.6 x lO" 2 r°^ d £)i r\ ~ 0.2 pc for the characteristic values of 3C 279. 

Alternatively, the acceleration process could become more efficient along the jet if the 
magnetic-field configuration gradually evolves into a more quasi-parallel one, on the same 
length scale of ~ 0.2 pc as estimated above. However, this scenario might be in conflict with 
the predominantly perpendicular magnetic-field orientation seen in the jets of 3C 279 on pc 
scales. 

The hard lag in the optical regime may also be indicative of a slow acceleration mecha- 
nism, with an acceleration time scale of the order of the observed B vs. R lag. This would 
imply an acceleration rate of 



lB-lR „ .. , rt _ 2 /obB \-l f D l 



1A ~ ^—^ ~ 6.8 x l0-\r^ d r ( jr- ) ' 2 s- 1 . (5) 
tbr \^gJ 

In this scenario, electrons could only be accelerated to at least 7#, if the acceleration rate 

of eq. [5] is larger than the absolute value of the radiative (synchrotron + Compton) cooling 

rate corresponding to eq. [3 This imposes an upper limit on the magnetic field: 

B < B acc ~ 0.42 D 1 / 3 (r^ d )- 2 / 3 (1 + A;)~ 2/3 G. (6) 

This can be combined with the estimate in eq. |4] to infer a limit on the magnetic-field 
equipartition parameter: 



e B < e B , acc ~ 5.8 x lO" 3 D?" « d )" 7 /3 (1 + fc) -r/3_ (7) 

Based on this equipartition parameter, one can use the magnetic-field estimate of eq. H] to 
estimate the total amount of co-moving energy contained in the emission region at any given 
time: 



K - U R 3 — - 2.5 x 10 49 D^T°% A (C s , d r 1/2 (1 + k ) e ^ ( 8 ) 

Assuming that the bulk of this energy is dissipated within the characteristic variability time 
scale, one can estimate the power in relativistic electrons in the jet: 



L jet ~ A ~ 4.5 x 10 43 < R s jd (Ctd)- 3/2 (1 + k) ergs s" 1 . (9) 

" var 
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Previous modeling works of the SEDs of FSRQs in general and 3C 279 in particular 
indicated characteristic magnetic field values of a few Gauss, in approximate equipartition 
with the ultrarelativistic electron population. The unusually low equipartition parameter in 
eq. [7] could therefore pose a problem for the slow-acceleration scenario. Note, however, the 
very strong dependence of e# on the Doppler factor (ex D 23 ^ 3 ). A Doppler factor D ~ 20 
could account for equipartition parameters of the order of one. Also, the energy requirements 
of eqs. [8] and [9] seem reasonable, and there does not appear to be a strict argument that 
would rule this scenario out. 

Another scenario one could think of would be based on a decreasing magnetic field along 
the blazar jet, leading to a gradually increasing cooling break in the underlying electron dis- 
tribution. This would require that the cooling time scale for electrons emitting synchrotron 
radiation in the optical regime would be equal to or longer than the escape time scale. Thus, 
the inequality in eq. [3] would be reversed. This would require unreasonably low magnetic 
fields. Furthermore, this scenario would be in conflict with the typically observed unbroken 
snap-shot power-law continuum spectra throughout the optical-IR range. Therefore, this 
idea may be ruled out. 

7. Summary 

We have presented the results of an optical-IR-radio monitoring campaign on the promi- 
nent blazar-type flat-spectrum radio quasar 3C 279 by the WEBT collaboration in January 
- April 2006, around Target of Opportunity X-ray and soft 7-ray observations with Chandra 
and INTEGRAL in mid- January 2006. Previously unpublished radio and optical data from 
several weeks leading up to the ToO trigger are also included. 

The source exhibited substantial variability of flux and spectral shape, in particular in 
the optical regime, with a characteristic time scale of a few days. The variability patterns 
throughout the optical BVRI bands were very closely correlated with each other, while there 
was no significant evidence for a correlation between the optical and radio variability. After 
the trigger flux level for the Chandra and INTEGRAL ToOs was reached on Jan. 5, 2006, the 
optical flux decayed smoothly by 1.1 mags, within 13 days, until the end of the time frame 
of the X-ray and 7-ray observations. The decay could be well described by an exponential 
decay with a decay time scale of = 12.8 d. The flux then recovered to approximately the 
pre-dip values in a much more erratic way, including a ~ 0.5 mag rise within ~ 1 d. 

A discrete correlation function analysis between different optical (BVRI) bands indicates 
a hard lag with a time delay increasing with increasing frequency, reaching ~ 3 d for the lag 
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of B behind R. This appears to be accompanied by a single indication of counterclockwise 
spectral hysteresis in a color- intensity diagram (B-R vs. R). Thus, spectral hardening during 
flares appears delayed with respect to a rising optical flux. There is no consistent overall 
trend of optical spectral hardness with source brightness. However, our data indicate that 
the source displays a rather uniform spectral slope of a ~ 1.9 at moderate flux levels 
(R > 14.5), while spectral variability seems common at high flux levels (R < 14.5). 

The occasional optical spectral hysteresis, in combination with the very steep IR-optical 
continuum spectral index of a ~ 1.5 - 2.0, may indicate a highly oblique magnetic field 
configuration near the base of the jet, leading to inefficient particle acceleration and a very 
steep electron injection spectrum. As the emission region propagates along the jet, a gradual 
hardening of the primarily injected ultrarelativistic electron distribution may be caused by 
the gradual build-up of hydromagnetic turbulence, which could lead to a gradually increasing 
contribution of second-order Fermi acceleration. This would imply a length scale of the build- 
up of hydromagnetic turbulence of Ar ~ 0.2 pc. 

An alternative explanation of the hard lag may be a slow acceleration mechanism by 
which relativistic electrons are accelerated on a time scale of several days. However, even 
though this model can plausibly explain the observed variability trends and overall luminosity 
of the source, it requires an unusually low magnetic field in the emitting region of B < 0.2 G, 
unless rather high Doppler factors of D 20 are assumed. Such a small magnetic field would 
be about an order of magnitude lower than inferred from previous analyses of simultaneous 
SEDs of 3C 279 and other flat-spectrum radio quasars with similar properties. 
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